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We develop the mean-field theory of a charge- density wave (CDW) state in magnetic field and 
study the properties of this state below the transition temperature. We show that the CDW state 
with shifted wave vector in high magnetic field (CDW^ phase) has at least double harmonic modu- 
lation on the most part of the phase diagram. In the perfect nesting case the single harmonic CDW 
state with shifted wave vector exists only in a very narrow region near the tricritical point where 
the fluctuations are very strong. We show that the transition from CDWo to CDW^ state below 
the critical temperature is accompanied by a jump of the CDW order parameter and of the wave 
vector rather than by their continuous increase. This implies a first order transition between these 
CDW states and explains the strong hysteresis accompanying this transition in many experiments. 
We examine how the phase diagram changes in the case of imperfect nesting. 

PACS numbers: 71.30.+h, 71.45.Lr, 74.70.Kn 



I. INTRODUCTION 

The properties of metals with a charge-density-wave 
(CDW) ground state attract great attention since the 
fifties (see, e.g., monographs [1,2]). In quasi-lD metals, 
the Fermi surface (FS) consists of two slightly warped 
sheets separated by 2kp and roughly possesses the nest- 
ing property that leads to the Peierls instability and fa- 
vors the formation of the charge- or spin-density-wave 
(SDW) state at low temperature. 

The mean-field description is known to be unable 
to describe strictly ID conductors, where the non- 
perturbative methods and exactly solvable models are of 
great importance. 3 However, in most materials manifest- 
ing CDW or SDW phenomena the nonzero value of the 
electron transfer integral between conducting chains and 
the 3D character of the electron-electron interactions and 
lattice elasticity reduce the deviations from the mean- 
field solution and also make most of the methods and 
exactly solvable models developed for strictly ID case 
inapplicable. 4 The effect of fluctuations in Q1D metals 
and their influence on the mean-field description of the 
Peierls transition in these metals has been considered in 
a number of papers (see e.g. [5,6] and references therein). 
It was shown that the interchain dispersion of electrons 
strongly damps the fluctuation and validate the mean- 
field description. In Q1D organic metals, 7 at which our 
present study is mainly aimed, the free electron disper- 
sion near the Fermi level is, approximately, given by 



e a {k) = hv F {\k x \ - k F ) - t±(k±) - aH, 



(1) 



where vf and kp are the Fermi velocity and Fermi mo- 
mentum in the chain (x) direction, H = (IbB is the Zee- 
man energy, [Lb is the Bohr magneton and B is the ex- 
ternal magnetic field. The perpendicular-to-chain term, 
is much greater than the energy scale of the 



CDW(SDW) transition temperature, T c0 . Only the "im- 
perfect nesting" part, t' b , of t±(k±) is of the order of T c o 
(see Eqs. (2), (3)). Hence, the criterion for the mean-field 
theory to be applicable, 5,6 t± 3> T c q, is reliably satisfied 
in most Q1D organic metals. 

The mean-field description of the CDW properties is, 
in many aspects, very similar to the BCS theory. The 
pairing of two electrons in superconductors is replaced 
in CDW by the pairing of an electron with the hole on 
the opposite sheet of the Fermi surface. The charge 
and spin coupling constants in CDW (see interaction 
Hamiltonian (5)) are analogous to coupling constants in 
spin-singlet and spin-triplet channels in superconductiv- 
ity. The CDW phase with shifted nesting vector is sim- 
ilar to the non-uniform superconducting phase (LOFF 
phase). 8 ' 9 However, there are several important differ- 
ences between these two many-particle effects. The first 
difference is that the formation of a gap in the electron 
spectrum in CDW leads to an insulating rather than su- 
perconducting state. This happens due to the pinning of 
the CDW condensate by crystal imperfections. Hence, 
the CDW state does not reveal a superfluid current. 1,2 
Other differences appear, e.g., in the excitation spec- 
trum. In particular, the lowest energy excitations in 
magnetic field or at imperfect nesting are not always the 
electron pairs as in BCS theory but may be the "soliton 
kinks". 10 - 14 

The theoretical investigation of the CDW/SDW prop- 
erties at a mean-field level comprises two main branches. 
The first focuses on the transition line from the metallic 
states to CDW/SDW state using susceptibility calcula- 
tions. This allows to include many additional factors into 
the theoretical model, such as different free electron dis- 
persion relations, spin and orbital effects of the external 
magnetic field, applied pressure etc. It helped to discover 
and to explain many beautiful effects such as the field- 
induced spin- and charge-density waves (FISDW), 18-21 



2 



the increase of the transition temperature due to the 
one-dimcnsionization of the electron spectrum in high 
magnetic field, 18 ' 23 the CDW state with shifted nesting 
vector 19 ' 24 ) etc. However, all these results cannot be 
continued below the transition temperature since they 
are based on the electron susceptibility calculations us- 
ing metallic-state electron Green's functions. In particu- 
lar, the calculation of [24] predicts an appearance of the 
CDW X phase with shifted wave vector as the magnetic 
field exceeds some critical value. This calculation gives 
the metal-CDW^; transition line, T C (H), and the depen- 
dence of the optimal shift q x (H) of CDW wave vector on 
magnetic field at this transition line. However, this calcu- 
lation does not extend below the transition temperature, 
and, hence, cannot be used to describe the properties and 
the phase diagram inside the CDW phase. For example, 
it neither describes the structure of the CDW X phase 
below T C {H) nor gives the CDWo-CDWa, transition line 
H C \(T) or the kind of this transition. The mean-field 
study of CDW in Rcf. [22] is applicable only in weak 
field when the CDW wave vector does not shift from its 
zero-field value. 

The second branch of investigation involves the soliton 
physics. It appeared first to describe the ground state of 
polyacetylene 10 ' 11 and then developed into a rather large 
activity (see Refs. [15,16] for a review). In particular, the 
soliton structure and the energy spectrum of CDW state 
in external magnetic field were considered 17 . However, all 
these results were only derived at zero temperature and 
perfect nesting and are not applicable at temperatures of 
the order of transition temperature, T ~ T c . The finite- 
temperature phase transition to soliton phase has been 
considered in Rcf. [14]. However, that phase diagram 
refers to zero field and finite shift of the chemical poten- 
tial from the value fiQ corresponding to the commensu- 
rate CDW. 14 The analysis of the CDW phase diagram in 
magnetic field has been performed in the case of perfect 
nesting and at the electron density close to half- filling. 25 
However, this approach is also unable to describe the in- 
commensurate case far from half-filling. In particular, 
it suggests the one-harmonic modulation of charge den- 
sity, which is usually unstable in high magnetic field (see 
below) . 

Further theoretical description of the CDW state in 
magnetic field became important last years because of 
the intensive experimental study of the CDW state in 
strongly anisotropic organic metals. 26-36 

In the present paper we study the CDW phase diagram 
and the properties of the CDW state at finite tempera- 
ture below the transition point in the mean-field approxi- 
mation. This study links the results of susceptibility cal- 
culations and the results from the soliton approach. We 
take into account the spin effect of the external magnetic 
field and consider the CDW phase with shifted nesting 
vector. The " antincsting" term in the electron dispersion 
is also taken into account. We show that the CDW X state 
with shifted wave vector, proposed in Ref. [24], has dou- 
ble harmonic modulation almost everywhere in the phase 
diagram (see Fig. 3,5). At perfect nesting, t' b = 0, the 



single harmonic CDW^ state exists only in a very nar- 
row region near the tricritical point. We show that the 
transition from CDWo to CDW X state below the criti- 
cal temperature is accompanied by a jump of the CDW 
wave vector rather than by its continuous increase. This 
implies a first order transition between the CDW states 
and explains the strong hysteresis observed at the kink 
transition. 28 - 29 - 32 ' 34 

Besides the Zeeman splitting, external magnetic field 
affects the orbital electron motion. First, an external 
magnetic field perpendicular to conducting layers leads to 
one-dimensionization of the electron spectrum in quasi- 
1D metals 23 that improves the nesting property of the 
Fermi surface. Second, due to the strong scattering by 
the nesting vector, electrons in quasi- ID metal in mag- 
netic field may form close orbits in momentum space 19 . 
These orbital effects increase the CDW or SDW transi- 
tion temperature and may lead to the field-induced spin- 
density waves 18 (for a review, see [7]). A similar effect 
for the charge density waves was also proposed 21 . The 
interplay between spin and orbital effects of the magnetic 
field is quite interesting. It leads, for example, to a se- 
ries of phase transitions between the CDW states with 
quantized nesting vector. 32 

In our present study we disregard the orbital effects 
of magnetic field. This limitation, however, is not very 
restrictive. If the magnetic field is parallel to the con- 
ducting layers, it produces only the Zeeman splitting 
and no orbital effect. The orbital quantization effects 
are, usually, more subtle than the spin effect. As an 
illustration, the quantized nesting phases can be ob- 
served only at very low temperature 32 because the ef- 
fects of nesting quantization are strongly damped at 
huj c = ehH/m*c < 2ttT, 2ttH/t, A. This damping is 
somewhat similar to that of the magnetic quantum oscil- 
lations. The one-dimcnsionalization of electron spectrum 
in magnetic field can be effectively taken into account in 
the present study by introducing the magnetic field de- 
pendence of the imperfect nesting transfer integral, t'(H), 
in the electron dispersion relation (2). A situation which 
is mathematically equivalent to the Zeeman splitting of 
magnetic field without its orbital effect arises when there 
are two slightly different Q1D chain system in the same 
compound with coupling between the chains of differ- 
ent type. The quantities U c + U s and U c — U s play the 
role of the electron coupling constant inside the chains 
and between the chains of different type, respectively. 
Such a system occurs in many compounds as TTF-type 
organic metals (see, e.g., the discussion on page 196 in 
[16]). Two slightly different types of chains occur also 
in a-(Per)2M(mnt)2- 35 ~ 37 In the systems with molecular 
chains of two types, however, the difference between op- 
timal CDW wave vectors on different chains is fixed by 
the molecular structure and does not vary as in the case 
of external magnetic field which gradually separates the 
Fermi surfaces with two different spin components. 
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II. THE MODEL AND THE MEAN-FIELD 
THEORY OF CDW IN MAGNETIC FIELD. 

1.. The model 

We consider a quasi- ID metal with dispersion (1) and 
t±(k±) given by the tight-binding model: 

t±(k±) = — 2tb cos(kyb) — 2t' b cos(2k y b) — 2t c cos(k z c z ), 

(2) 

where b and c z are the lattice constants in y- and z- 
directions respectively. The dispersion along the z-axis 
is assumed to be much weaker than the dispersion along 
y-direction. Therefore, we omit the second harmonic oc 
cos(2fc z c z ) in the dispersion relation (2). Since the terms 
2t c cos(k z c z ) and 2tbCos(k y b) do not violate the perfect 
nesting condition 

e(p + Q) = -e(p), (3) 

they do not influence the physics discussed below unless 
the nesting vector become shifted in y-z plane. We do 
not consider such a shift in the present study. 
The electron Hamiltonian is 

H = H + i?int , 

with the free-electron part 

i?0 = 5>a(fc)4(fc)a CT (fc) (4) 
ka 

and the interaction part 

Him = I V aa ,(Q)ai(k + Q)a a (k) 

kk' Qaa' 

xal(k' -Q)ay{k'). (5) 

This Hamiltonian does not include the orbital effect of 
magnetic field on quasi- ID electron spectrum which may 
be important in the case of substantially imperfect nest- 
ing and strong magnetic field perpendicular to the easy- 
conducting plane (see introduction). Later on we will be 
interested mainly in the interaction at the wave vector Q 
close to the so-called nesting vector Q , which is usually 
chosen as 

Q = {±2k F ,Tr/b,n/c). (6) 

The deviations Q — Q that will be considered below are 
of the order of max{7J, t' b }/hvp <C k Fl so that for such 
small deviations the interaction function 

V a u 1 (Q) « Vw'iQ = Q ) = U C - U s aa' (7) 

The coupling constants U c and U s are the same as in 
Ref. [24]. Subscripts c and s distinguish charge and 
spin coupling constants. The Hamiltonian (5) can be 
obtained from the extended Hubbard model when the e- 
e scattering by the momenta close to the nesting wave 
vector Qo is only taken in to account. 



2.. Mean field approach 

First, we note that the terms with Q — in (5) only 
renormalizes the chemical potential and will not be omit- 
ted in subsequent calculations (all sums over Q do not in- 
clude Q = 0. We introduce the thermodynamic Green's 
function 

g a (k', k,T- t>) = (T T {4(fe' ; r'K(fc, r)}>, (8) 

where the operators are taken in the Heisenberg repre- 
sentation, and the average 

D Q ^J2^( k -Q^^-°)- s (Q^) n - ( 9 ) 

k 

After defining 

Aq<t = Y,Vc - U s aa')D Qa , (10) 

a' 

one obtains in the mean-field approximation 
Hi nt = a t( k + Q)a*(k)A Qa - g 2 £) -Q- A Q- 

Qk<7 Qcr 

(ii) 

Hermiticy of the Hamiltonian requires A_q ct = Aq ct . ft 
is now straightforward to write down the equations of 
motion which in the frequency representation take the 
form 

[iu - e a (k)]g a (k', k,u))~y~] A Qa g a (k', k-Q,w) 

Q 

= 6{k',k). (12) 



3.. The Cosine Phase 

Now we consider the solution with Afc CT ^ only for 
k = ±Q, where Q — 2k F e x + (ir/b)e y + q with \q\ -C kp. 
If we neglect the scattering into the states with \k x \ > 
2kp, the equations (12) decouple: for k x > one has 

V -A_ Q(T iw-e a (k-Q)J v ; 

where 

r = f 9a{k,k,uj) g a (k - Q,k,w) \ , . 
U -\g a (k,k-Q,Lj) g a (k-Q,k-Q,u)J 

and / is the 2x2 identity matrix. The equation for 
k x < may be obtained via a substitution Q — > —Q. 
Introducing the notations 

et(k\k) = £ ° {kl) ± £ ° {k) (15) 

and using AQ a A-Q a = |Aq ct | 2 from (13) one has 
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g<r(k - Q,k,w) = 
The consistency equation therefore is 



[lu + ie a (k)][w + is a {k - Q)] + A Qct | 2 

Aq^ 

" [w + fe - Q)]2 + [ e -(fe, fe - Q)]2 + |A Qa p ' 



(16) 



(t/ c - £W)A Q(T 



+ 2e+(fe, fe - Q)] 2 + [e^ (fc, k - Q)f + \A Qa \ 



(17) 



where uu takes the values nT(2n +1), n € Z. Using the identity 



i 



(w — ia) 2 + a/ 



2T 



th 



a + a 
2T 



this may be rewritten as 
where 



4e 



fc' ^/|A Q(T | 2 + Mfc,fc-Q)] 2 ^ 2T 2T 



E ai± {k) = e+(k, k-Q)± y |A Qff | 2 + [e^ (fc, fe - Q)] 2 
is the energy spectrum in the cosine phase. From (2) it follows that for the right-moving electrons 



£+(fc, k — Q) = ^ V ^ x -f ■>!,. sin sin ; /,-„/> - ----- 



(kyb ~ — cos <fe^ cos (2fc y 6 — <lyb) — crH, 



e a (k, k — Q) = hvp(k x — kp — q x /2) — 2% cos cos ^k y b — + 2t / b sinq y bsm (2k y b — q y b) . 



(18) 

(19) 

(20) 
(21) 



It is, however, more convenient to perform the momentum summation in Eq. (17) first. Extending the limits of the 
summation over k x to infinity, one has 



nv F \Uc\T \ - 
Aq ct < ^ 



(1 + vaa')A Qo 



[io + ie+{k,k-Q)f + \A Qo 



(22) 



I 

where the branch of the square root with positive real (the superscript "c" stands for "cosine") with 
part is implied, the angular brackets stand for the averag- 
ing over all values of k y , and we introduced the standard Ty . m _ lTT lrT1 ^ / signw 
notation 



v = -U s /U c 
for the coupling ratio and 

v F = 



nhvF 



(23) 



(24) 



KW = ku f \U c \tY, 



uj + iet(k, k Q) , fe< 

ttvf\U c \T <sr-^ I signw 

(oj + iei{k,k-Q)f I k 



(26) 



(27) 



for the density of states on the Fermi level per one spin 
component. 

Expansion to the third order in Aq ct yields 



The second-order transition line corresponds to 

= 0. (28) 



dct I V< 1} ^> - 1 



Aq,j' = 1 + wo — \Aq„\ 2 K^ c \ Aq ct (25) The left hand side of this equation is a function of is, T, H 

2 V " a J anc l q. The normal to CDW phase transition occurs at 
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some optimal value of q which corresponds to the min- 
imum of the free energy of the CDW state, that at the 
transition point is equivalent to the maximum of T c or to 
the minimum of l.h.s. of Eq. (28). 

Unrestricted summation over k x introduced divergence 
into the summation over u> in which may be renor- 

malized by introducing the zero-field and zero-i^ transi- 
tion temperature T c q. To obtain a closed expression for 
the latter, we note that for H = one has = 
and the matrix in (28) has eigenvalues — 1 and 



vK { „ ] - 1. Therefore, at H = K$\t = T c (t' b ),H = 
0) = 1. From (26) at t' b — after imposing a cutoff at 
Ld ~ E F one obtains if (1) = v F \U c \ \n(2-fE F /nT), and 

T c0 = ^cx P (-4-), (29) 

7T V V F \U C \) 

in agreement with Ref. [39] and with the usual BCS 
expression. We may now write 



l + u F \U c 



In- 



T, 



c0 



signw 



w + iej(k,k - Q) I k w 



signw 



(30) 



The kernels and K^ ,c can be further simplified in terms of the digamma function ip(x) = 4- lnr(x): 



,(3)c 



K<P = \ + vf\U c 



1 Tc0 I / 
In — + V 



Rc?/> 



1 ie+(k,k-Q) 

2 2ttT 



vf\Uc\ 
16tt 2 T 2 



Rei/*" 



„ (I , is+(k,k-Q) 



+ 



2irT 



(31) 
(32) 



To actually find the fields Aq ct just below the transition point to the leading order, one has to make the substitution 



-a 

(28), 



Ky — > K)j } — \/S.Q a \ 2 K^ c in (28). The ratio Aq_ : Aq + is real and at the transition point is the column ratio in 



: A 



Q+ 



I+2^« i 



(33) 



Introducing the (real) ratio a = |Aq_| 2 : |Aq + | 2 we get 

a (i-)/2 (vK^kW - ^(K™ + K&) + 1 



(34) 



The superscript "c" stands for the cosine phase. The 
CDW cosine phase with q x ^ was analyzed at the tran- 
sition line by means of the susceptibility calculation and 
called the CDW^ phase, while the phase with q x = was 
called CDW 24 . 



4.. The Double Cosine Phase 

Now we consider the solution with Ak a ^ for k = 
±Qo ± <7, where Qo = 2k F e x + (jr/b)e y and q <C fc.F- 
Strictly speaking, there exist no self-consistent solutions 
with only four harmonics present or all others damped 
by a factor of AQ a /E F . However, the modulations given 
above have equal second-order transition temperatures, 
and immediately below the transition point other har- 
monics are damped by a factor of Aq,j/T or Aq^/H. 

To obtain the leading-order expressions for the fields 
AQ a in the vicinity of the transition line, we rewrite the 



equations (12) in the "matrix" form 



GGq 1 -GF = I, 



(35) 



where 



G(k',k) = g a (k',k), G (k',k) = 

F{k',k) = Y J ^Q^(k',k-Q). 

Q 



6{k',k) 
iu - e a {k) ' 



(36) 



The solution is 

G = Go + G FGo + G FG FGo + ... (37) 

Omitting the contribution from the virtual states with 
the momentum \k x \ > 2k F , we obtain the consistency 
equation to the third order in Aq ct : 
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A Qo+q . a , =tJ2(Uc-U s <tct') 

kujcr 

+ 



A 



Qa+q,(7 



(iuj - s a (k))(id> - e a (k - Q Q - q)) 



^Qo+qM^Qa+qM 



(iu> - s a (k))(iu - e a (k Q Q q)){iu - e a {k)){iuj - e a (k Q Q q)) 



+ 



+ 1- 



A 



(38) 



Qo+q^^Qo-qM 



+ 



(iu - s a (k))(iu> - e a (k - Q - q))(iu> - e a {k - 2q))(iu> - e a (k Q q)) 

^Qa+gA^Qo-qA 2 

(ioj - e a (k))(iu - e a (k - Q + q))(iw - £ a (k))(iu - s a (k Q a q)) 



+ 



The equation for the other two harmonics may be obtained via a substitution q — > —q. Performing the integration 

1 + van' 



over k x , we arrive at 



a 

(the superscript "d" stands for "double cosine" ) with and K^ c given by (31) and (32), and 



(39) 



TT VF \U C \T^2 

to 

V F \Uc\ I 

A-Khv F q x T \ 



(u> — iaH) sign a; 



[(w + ie+{k, k - Q))(w - iet a (k, k - Q))] s 



Im ip' 



,(\ , ie+(k,k-Q] 



2irT 



Imip 



,/l ie± a (k,k-Q) 
2ttT 



(40) 



The function K„ does not depend on a and in what 
follows we omit this subscript. 

At only a longitudinal shift of the CDW wave vector 
(q y = 0) the fc y -dependence of e+ a (k,k — Q) is sym- 
metric, and the functions K^ c and K^ d (being 
dependent only on the shift wave vector q x ) possess the 
symmetry: 

K^{q x )=K { *l C {-q x l (41) 
K^ d {-q x )=K^ d {q x ). 

This symmetry follows directly from Eqs. (31), (32) 
and (40) and the properties of digamma function: 
Reip^ (a + ib) = ReV (n) (a-«&), lm^(a + ib) = 



— Imip( n ^ (a — ib). This symmetry is not a consequence 
of the expansion in powers of A CT : the consistency equa- 
tion (22) also does not change under the transformation 

Qx > 9x) & > (42) 

The ratio |AQ 0+qj(T | 2 : \AQ ()+q _- a \ 2 near the tran- 
sition point is still given by (33) accounting for the 
property (41). Hence, |A Qo _q ;(T | 2 : |A Qo _ q _ CT | 2 = 
lAQo+q^l 2 : \A Qa+qtCJ \ 2 = a. As in derivation of (34) 
one may rewrite the system of equations (39) on A as 
two equations (28) with replacement K^~\Qq ± q x ) — > 
K { H \A Qo±q ^\ 2 K^ c \A QoTq ^K^. 

Defining the ratio (3 = |AQ _ q ,_| 2 : \A Qa+q , + \ 2 we can 
extend the derivation of (34) to get 



^Qo+q,+ \ 



and from the Qo ~ Q part of (39) 

,(i)^(i) _ 



(43) 



|A Qo - 



(44) 
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Dividing (43) by (44) we obtain a linear equation on (3: 



(45) 



which gives (3 — 1. The values /3 = or /3 = oo in 
(43) and (44) correspond to cosine phase. Hence, in the 
double cosine phase the symmetry (42) is not broken, 
while the transition to cosine phase breaks this symmetry. 



5.. Free energy of cosine and double-cosine phases 

One can easily write down the expressions for the free 
energies of cosine and double-cosine phases valid to the 
second order in the energy gap |A CT |. 

From (11) we have for the free energy 



The quantity + aK^' depends on q but is always 
positive near the metal-CDW transition line T C (H) where 
this transition is of the second kind. 

Since CDW C and CDW2 C phases have the same transi- 
tion temperature, the only way to determine which phase 
takes place is to compare their free energies, that near the 
transition line T C (H) are given by the formulas (48) and 
(49). The double cosine phase wins if the ratio 



-(I) 



-FCDW - F n = - ^2 D-QctAqv. 

To the second order in Aq ct this rewrites as 

ttvfT ^ I AqctA-Qct signw 



(46) 



-Fcdw — F n = — 



E 



Qu* ^ + ie7(k,k-Q)/ kv 



The phase with the most negative r.h.s of (46) wins; pos- 
itive values of the r.h.s. correspond to first-order transi- 
tions. From Eq. (47) using Eq. (41) we have 



P TP - + +Q:iV - I A C |2 

c n ~ 2\Uc~\ Q °+ q ~ +l 



(48) 



and 



"2c 



_ + -raj\_ |a2c 2 
™ - \u I I^Qo+qi,+ l • 



r F 



"2c 



F n ) 



(F c - F n ) 



> 1, 



(50) 



where the values of the functions F% C (T, H, q) and 
F C (T, H, q) must be taken at the optimal value of the 
wave vector q, that should be found by minimization of 
these free energy functions at each point of T, H phase 
diagram and for each of the two phases (CDW C and 
CDW2c) separately. 

Below the transition temperature the optimal shift vec- 
tors q could be different for the cosine and double cosine 
phases (they minimize different free energy functions F c 
and F2 C )- However, at the transition temperature T C {H) 
the CDW C and CDW2 C phases have the same optimal 
value of q which is determined by the minimum of the 
left-hand side of Eq. (28). Hence, on the transition line 
T C (H) Eq. (50) simplifies to 



r F = 



o|A2c 12 

I A c I 2 

I^Qo+q x ,+ 



> 1, 



(51) 



(49) or after substitution of (34) and (43) 



r F 



> 1. 



(52) 



These formulas will be used later to determine the phase 
diagram. 



III. THE PHASE DIAGRAM 

In this section we only consider the longitudinal mod- 
ulation (q y = 0) of the CDW wave vector. The CDW 
phase with q y ^ (CDW y phase) may appear for cer- 
tain dispersion functions t±(k±). However, for the tight- 
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binding model (2) the CDW y phase is not expected to 
take place near the transition temperature. Near the crit- 
ical pressure (when t' b « t' b * and the SDW-metal phase 
transition takes place) in zero field the SDW2 phase was 
predicted using the susceptibility calculation from the 
normal state. 40 However, this phase was predicted only 
near T = 0. Moreover, the absolute instability line cal- 
culated in [40] does not give the actual transition line at 
low temperature since this transition takes place as the 
first-order phase transition. 

To be more sure that we can disregard appearance of 
the CDWj, phase in our consideration near the transition 
temperature we performed the calculation of the optimal 
shift vector q at the transition line T C (H) in a large range 
of parameters. We swept the following three-dimensional 
range of parameters: —0.3 < v < 0.9,0 < t' b /t' b * < 1 and 
< hbH < 2T c q and have not found that the shift of Q y 
leads to higher T c anywhere inside this range. 

This does not contradict the prediction of the CDWy 24 , 
since according to this paper the CDW y phase appears 
at v < essentially due to the orbital effect of the mag- 
netic field (see Fig. 7d of Ref. [24] ) . Below we will only 
consider the range < v < 1 which is close to the ex- 
perimental situation in organic metals and disregard the 
appearance of CDW y phase. 



1.. Transition line T C (H) and the 
normal-CDWo-CDWt tricritical point 

The normal-to-CDW phase transition line is given by 
Eq. (28) irrespective of to which CDW phase this tran- 
sition occurs. The shift of the wave vector q x which en- 
ters Eq. (28) corresponds to the maximum value of T c 
at given magnetic field H . The behavior of T C (H) and 
q x (H) at the transition temperature has been analyzed 
in Ref. [24]. Diagonalization of the M^^M\ part of the 
susceptibility matrix (6) in Ref. [24] corresponds to the 
diagonalization of the matrix in Eq. (28). On the transi- 
tion line one of the eigenvalues of these matrices is zero. 
In the case of perfect nesting equation (28) gives the same 
result for T C (H) as Eq. (17) of Ref. [24] (sec fig. 7(a) of 
[24]). 

For q x = (CDW phase) K ( + ] = K W and Eq. (28) 

simplifies to = 1, that gives for the transition tem- 
perature T c = T C (H) the well-known equation 



where 



In 



T c (H,f b ) 



i_z £ +(fc,fc-Qo) 

2 27rT c 



(53) 



e+(fe, k - Qo) = -2t' b coa(2k y b) - <jH. (54) 



In the Eq. (53) one can take either of the values a = ±1 
since this equation does not depend on the sign of a. The 
transition line T C (H) from normal to CDWo phase does 
not depend on the value of v. 

The tricritical point (where the CDWo and CDW^ 
phases have the same transition temperature) is given by 
equation d 2 D(v,T,H,q x )/dq x = 0, where D(v,T,H,q x ) 
is the left-hand side of Eq. (28). At this point q x = 0, 
K { + ] = K a) = 1 and d n K ( ^ ) /dq r x l = {-i) n d n K W /dq™, 
and the equation for the tricritical point becomes 



02^(1) 



2v 



(dq x ) 2 1 - v \ dq x 
Substituting (31) and (32) this equation rewrites as 



-( ReV> 



M / 1 ie+(k,k-Q Q ) 



2ttT c 



where 



, I 2 27TT, 



V = 2vv F \U c \/(l-v). 




(56) 



(57) 



Together with Eq. (53) this equation allows to find the 
tricritical point 



2.. T C (H) transition line at perfect nesting 

The case of perfect nesting (t' b = 0) is equivalent to the 
strictly ID case in many mathematical aspects. However, 
the 3D nature of the compound (its energy spectrum, 
e-e interaction and lattice elasticity) preserves, and the 
mean- field approach still works. The analysis in this sub- 
section (at t' b = 0) is simple and performed analytically. 
It helps to understand some qualitative features of the 
CDW phase diagram in magnetic field. 

At perfect nesting and purely longitudinal modulation 
of CDW (q y = 0) the expressions for K^ d and 

K^ c simplify to 
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K% = l + v F \U c 



In- 



T, 



cO 

T 



K, 



(3)d 



vf\U c 



4Trhvpq x T 



1 2 2ttT J 



Reip 
ih„ \ 



iha_\ 
2nT J 



, i 1 ih- a 

+ Imib' — I - 

1 2 2fT 



(58) 
(59) 



and 



where 



ff ° 16tT 2 T 2 



W^ + ^J. (60) 



(61) 



At g x = 0, Eq. (59) even more simplifies: 

"(l-^\= 2K?H (62) 



K (3)d_ VF\U Ui R ,, 



One can now derive a simple formula for the transition 
line T C {H) in the high field limit H » T C {H). As H ^ 
oo, in (61) goes to zero for one spin component and 
to —2aH for the other. Using the limit expansion of the 
digamma function, Re -0(1/2 + ix) = Inx + 0(l/x), x — > 
oo, from (58) we have at H/nT c (H) > 1 



*l + v F \U c \ln(T c0 /T) 
K ( y » 1 + ^ F |*7 C | In (7rT c0 /4 7 i/) . 

Equation (28) rewrites as 



(63) 



(l + z/)if a) /2-l 
i/Jfi 1} -(l + i/)/2' 

that in the limit H/ttT c (H) » 1 becomes 

ln(4 7 ffAT c0 ) 



In(Tco/T) 



^111(47^/7^00) + 1' 
At ?7ln(47i?/7rT c0 ) > 1 this simplifies to 

T C (H ^00) =T c0 exp(-l/r?). 



(64) 



(65) 



In the intermediate interval ttT c (H) -C 47 <C 
7rT c0 exp(l/77) we get 



T c (i?) w 7rT c 2 /(4 7j ff). 



(66) 



Hence, at perfect nesting in the limit H — > 00 the transi- 
tion temperature tends to a finite value which is natural 
since it corresponds to the CDW instability for only one 
spin component. 

However, the behavior of T C (H) given by formulas (64), 
(65) and (66) strongly changes at finite " antincsting" 
term t' b (see below). 



3.. Metal-CDW transition lines and tricritical 
points at finite t' b 

Using the Eq. (28) and Eq. (56) we calculated the 
transition lines T C (H) and the tricritical points at differ- 
ent values of t' b and v. The results are shown in Fig. 1 . 
From this figure we see that the transition line T C (H) de- 
pends strongly on the both parameters t' b and 1/, and one 
cannot determine the value of t' b from the T C {H) transi- 
tion line if one does not know the value of v. Hence, to 
determine experimentally what values take t' b and v in a 
particular compound, one has to perform two indepen- 
dent tests. For example, one can study the pressure de- 
pendence of the T C (H) lines at two different tilt angles of 
magnetic field: the first measurement when the magnetic 
field is parallel to the conducting x-y plane and the for- 
mulas derived above in neglect of orbital effects are valid, 
the second measurement in the very strong magnetic field 
perpendicular to the conducting x-y plane, when the or- 
bital effects are so strong that the one-dimensionization 
of the electron dispersion takes place and the formulas 
for perfect nesting (see sec. 3.2) become approximately 
valid. One can also determine the value of t' b indepen- 
dently from the transport measurement and then com- 
pare it to the critical value of t' b * = Ao/2 « 0.88T c o at 
which the CDW state is damped without magnetic field. 

There is one interesting common feature on all dia- 
grams in Figs. 1. At each value of v there is a critical 
value of t' b < t b above which the CDW X phase disap- 
pears, i.e. as magnetic field increases the transition from 
CDWo to metal state instead of the CDW X phase takes 
place up to the lowest temperature. This critical value 
of t' b may shift due to the "one-dimensionization" effect 
of magnetic field. 

We have also checked if the region of cosine phase 
increases considerably after finite " antinesting" term is 
taken into account. The results of this study at the tran- 
sition line are given in Fig. 2. These results indicate that 
the region of cosine phase does not increase considerably 
at finite t' b for the tight-binding dispersion. 



4.. Cosine and double-cosine phases at perfect 
nesting 

To determine which of the two phases (CDW C or 
CDW2c) wins we have to compare their free energies 
given in section 2.5. 

On the transition line this reduces to the evaluation of 
the ratio (52). At the tricritical point (when q x = 0) this 
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Magnetic field (Hn B /T co ) 

FIG. 1: Transition lines T C (H) at different values of t' b at three 
different values of the coupling constant ratio v. The values of 
t' b is given on the figures in units of T c o; for the tight-binding 
dispersion (2) the critical value of t' b at which CDW disappears 
without magnetic field is t' b * = Ao/2 « 0.88T cO . The triangles 
show the tricritical points Normal-CDWo-CDW^. 



1.3-, 




I 



0.0 0.2 0.4 0.6 0.8 1.0 

v=U s /U c 

FIG. 2: The magnetic field values H c \(v) and H 02(f) at two 
tricritical points Normal-CDWo-CDWc and Normal-CDW c - 
CDW2C correspondingly as function of coupling constant ratio 
v — Us/Uc for two different values of t' b . The solid and dot 
lines denote H C 2(v) and H c \(v) at perfect nesting (t' b = 0) 
while the dash and dash-dot lines denote H C 2(v) and H c i(i/) 
at finite value of the "antinesting" term (t' b = 0.5T c o). The 
cosine phase exists in the areas between H c \(v) and H C 2(f) 
lines at the same value of t' b . 

ratio can be easily evaluated using Eq. (62): 

r F tr ip ie = 2/3 < 1. (67) 

Hence, near the tricritical point the cosine phase wins. 
In the limit of high field (H » T c ) one has hvpq x — > 

2H and h„ -> (1 - a)H. In this limit K { + )c w 

-(^ F |[/ c |/167r 2 T 2 )RcV"(l/2) > Xi 3)c , X ( 3 ) d , and from 
(52) we get 

r F -> 2 at H/T c -» 00. (68) 

Hence, at high magnetic field the double cosine phase 
wins. These two simple estimates suggest that the co- 
sine and double-cosine phases both appear on the phase 
diagram. 

The boundary between CDW C and CDW2 C phases 
on the transition line T C (H) is given by the equation 
Tf(T,H,v) = 1, which rewrites as 




This equation is valid also for nonzero t' b and together 
with Eq. (28) allows to determine the second tricritical 
point (H C 2(v),T c (H C 2)) where the normal, CDW C and 
CDW2C phases meet. 

In Fig. 2 we plot two triple-point values H c \(v) and 
H ci(v) of magnetic field which are given by equations 
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CDW 



normal phase 




0.55 



n — ' — r 

0.90 0.92 0.94 0.96 0.98 1.00 1.02 1.04 1.06 1.08 1.10 

H/Tc 

FIG. 3: Phase diagram in H — T coordinates at v = 0.7 at 
perfect nesting near the tricritical point. The CDW C region 
disappears rapidly as temperature decreases below T c . 



(56) and (69) correspondingly (together with Eq. (28)) at 
two different values of t' b = and t' b = 0.5T c0 . H cl {v) and 
Hciiy) determine the tricritical points Normal-CDWo- 
CDW C and Normal-CDW c -CDW 2c - From this figure we 
see that the range of the cosine phase is very narrow (the 
difference H c2 — H c i does not exceed 5% of H cl ) and 
strongly depends on v. 

At perfect nesting two tricritical points H c \ and H C 2 
coincide at v — 0. At this "double-triple" point the de- 
nominator in formulas (34) and (43) is zero. This means, 
that one should expand up to the higher orders in |A 2 | 
than fourth and that at this point the critical fluctuations 
are very strong. 

From Fig. 2 we see that (i) on the phase diagram of 
CDW in magnetic field there are at least 3 different CDW 
states: CDWo, CDW C and CDW2 C , (ii) at the transition 
line T C (H) the region of the CDW C phase is quite nar- 
row for various values of t' b and is usually sandwiched 
between the CDWo and CDW2 C phases as the magnetic 
field increases. 

To analyze how this picture evolves below the tran- 
sition line T C (H) we performed a numerical calculation 
of the free energies of CDW C and CDW2 C phases using 
formulas (48), (49), (43), (34). These formulas are valid if 
A <C ttT, H, that covers a narrow region below T C (H). 

The computation performed in the case of perfect nest- 
ing shows that the region of the CDW C phase becomes 
even more narrow as the temperature decreases and dis- 
appears at T w 0.95T C (H) (see Fig. 3). Hence, the cosine 
phase at perfect nesting exists only in a very small region 
of the phase diagram where it is strongly smeared by the 
critical fluctuations. 

The similar phase diagram appears in surface super- 
conductors in parallel magnetic field. 41 The cosine and 
double-cosine CDW phases correspond to the helical and 
cosine (stripe) superconducting phases respectively, and 



0.50 




0.45 



O 
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8 40 
> 

l 0.35 

Q 

O 0.30 



0.25 



Tricritical point 
normal-CDW -CDW/ 

C 2C 



0.60 0.62 0.64 0.66 0.68 0.70 0.72 

Temperature T/T 



CO 



FIG. 4: The shifts q x of wave vectors of the phases CDW2 C 
(solid line) and CDW C (dash line) along the CDW c -CDW 2c 
transition line (at 0.713 > T/T c0 > 0.66) or CDW -CDW 2c 
transition line (at T/T c o < 0.66) as function of temperature. 
These shifts coincide only at the normal-CDW c -CDW2c tri- 
critical point at T/T c o = 0.713. The difference between the 
solid and dash lines gives the jump of the CDW wave vector 
on the CDW c -CDW2c transition line that makes this tran- 
sition of the first kind and leads to hysteresis. The dash-dot 
line stands for the normal-CDW transition temperature, and 
the dash-double-dot line shows the temperature of CDWo- 
CDW c -CDW 2c tricritical point. 



the total CDWz phase corresponds to the non-uniform 
superconducting LOFF state. 8,9 However, our situation 
differs from that in Rcf. [41], where the Rashba term in 
the electron dispersion relation plays an important role 
in obtaining the phase diagram similar to that on Fig. 3. 
The number of harmonics and coupling constant in our 
case is twice larger than that in Ref. [41]. However, the 
analogy between CDW and surface superconductors in 
magnetic field is rather deep (see the discussion section) . 

On the transition line CDW C -CDW2 C the energy gaps 
and the optimal shifts q x of the CDW wave vectors of 
these two phases become substantially different as the 
temperature decreases below T c (see Figs. 4). This 
means that the transition from CDW C to CDW2 C is of 
the first order. The transition line between CDWo and 
CDW2C is also of the first order, and the CDW2 C nesting 
vector differs from Q — 2kp already on the transition 
line (see Fig. 4). This fact explains the huge hysteresis 
in magnetization and magnetoresistance observed in a- 
(BEDT-TTF) 2 KHg(SCN) 4 at the kink transition. 29,30 - 32 
The first-order phase transition driven by magnetic field 
is usually accompanied by hysteresis. The jump of the 
nesting vector on the CDW0-CDW2C transition line itself 
leads to huge hysteresis. The CDW and its wave vector 
are pinned by the impurities and crystal imperfections. 1 
As magnetic field increases through the kink transition 
point H = H C (T), the jump of the CDW wave vector 
forces the CDW condensate to move. Due to the pinning 
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this motion is hysteretic since at each time moment the 
CDW condensate finds some local minimum of the impu- 
rity potential. The larger jump of the CDW wave vector 
at the kink transition the greater the hysteresis is. This 
jump decreases with increasing temperature and comes to 
zero on the transition line T C (H). The hysteresis shows a 
similar temperature dependence. The hysteresis reduces 
with heating also due to the thermal activation processes 
which reduce the pinning. 

IV. DISCUSSION 

The electron charge and spin densities and the lat- 
tice distortion differ considerably in the cosine CDW C 
and double-cosine CDW2 C phases that can help to dis- 
tinguish these two states experimentally. In the double- 
cosine phase the charge modulation is the sum of two 
cosine distortions that leads to the beats of the charge 
density wave: 

p${x) = u F {A + + A_){cos[(Q + q x )x + fa] 
+ cos[(Q - q x )x + fa]} 
= 2i/ F (A+ + A_) cos[Q a: + {fa + fa)/2] 

xcos[q x x + (fa-fa)/2}. (70) 

The phase shifts fa and fa may depend on coordinate 
because of the pinning of CDW by impurities. Usually, 
A+(<2) and A_(Q) in CDW2 C state differ substantially, 
and the charge density wave in magnetic field is accompa- 
nied by a spin density wave. The spin-density modulation 
in the double-cosine phase is given by 

pf(x) = v F (A+ - A_){cos[(Q + q x )x + fa] 
- cos[(Qo - q x )x + fa]} 
= -2u F {A + - A_) sin[Q a; + {fa + fa)/2] 

xsm[q x x+{fa-fa)/2}. (71) 

In the cosine phase both charge and spin densities have 
one-cosine modulations, 

Pc(x) = v F (A + + A_) cos[(Q + Qx)x + fa] 

p c s {x) = v F {A+ - A_) cos[(Q + q x )x + fa]. (72) 

Thus, the charge modulations in the cosine and double- 
cosine phases can be experimentally distinguished by X- 
ray or Raman scattering, and the spin modulation can 
be detected by muon or neutron scattering experiments. 

The energy spectrum in the cosine and double-cosine 
phases also differ strongly. In the cosine phase the en- 
ergy spectrum (19) is asymmetric with respect to the 
spin components since the energy gaps A a differ for two 
spin components a. This means that under external elec- 
tric field the spin current is produced in addition to the 
charge current since the charge is transferred by elec- 
trons with predominantly one spin component. The de- 
gree of current polarization depends on the shift of CDW 
wave vector, and, hence, can be controlled by the exter- 
nal magnetic field. This property of the cosine phase may 



find applications in spintronics. The doublc-cosinc phase 
is symmetric in spin components, and its energy spec- 
trum has at least two gaps for each spin component. The 
symmetry (42) preserves in the double-cosine phase while 
in the cosine phase it is spontaneously broken. Thus, 
the one-cosine and double-cosine CDW states differ sub- 
stantially in their thermodynamic, transport and optical 
properties. 

At low temperature T <C T c the expansion in (37) is 
not applicable and the additional harmonics Q + (2n + 
l)q x with integer n appear in the double-cosine solu- 
tion of the consistency equation. These harmonics make 
the charge modulation in the CDW2 C phase at low tem- 
perature essentially nonsinusoidal in space and, possi- 
bly, containing soliton walls. At low temperature the 
quasi-particles differs substantially from the pair excita- 
tions with activation energy 2 A suggested by the simple 
mean-field description. Even at finite inter-chain electron 
dispersion the variational methods using the soliton solu- 
tions of the consistency equations show 16 that the lowest- 
energy excitations do not resemble the one-electron-type 
quasi-particles but involve many electrons and are accom- 
panied by a " soliton" kinks in the coordinate dependence 
of the order-parameter A a (x). These excitations in the 
incommensurate CDW have spin 1/2, zero charge and the 
energy 2A/ir for a ID chain without magnetic field. In 
external magnetic field their activation may become ener- 
getically favorable at H > 2A /n 17 that corresponds to 
the phase-transition from CDWo to the CDW^ phase. 
Unfortunately, there is no a complete solution of this 
problem at nonzero temperature at present time. The 
calculation at the electron density close to half-filling 25 
does not describe properly the incommensurate case far 
from the half-filling. In particular, it considers only the 
one-cosine modulation of the charge density. As we have 
seen, it is usually energetically less favorable than the 
double-cosine phase in high magnetic field. Therefore, 
we do not present an exact phase diagram in the whole 
temperature and magnetic field range. 

Qualitatively, for the tight-binding dispersion one may 
expect the following picture. The cosine phase does not 
appear on the phase diagram besides the small region 
near the tricritical point as shown on Fig. 3. One can 
easily show that in the case of perfect nesting at T — > 
the cosine phase is always unstable toward the forma- 
tion of the double-cosine modulation. However, an ob- 
servation of the cosine phase would be very interesting 
because of its spin asymmetry in thermodynamic and 
transport properties. For this reason some further calcu- 
lation of the phase diagram for various electron disper- 
sion relations would be interesting. The transition from 
CDWo to CDW2 C phase at low temperature comes about 
as the appearance of soliton kinks (or soliton walls). 17 
The soliton phase described in Rcf. [17] has two for- 
bidden bands (energy gaps) which makes it similar to 
the double-cosine phase. With increase of magnetic field 
the density of kinks increases, and the order parameter 
comes continuously to the sinusoidal modulation close to 
that in the double-cosine phase. At perfect nesting the 
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FIG. 5: Phase diagram in H~T coordinates at perfect nesting 
in the whole essential region of CDW at v = 0.7. The solid 
lines have been calculated from Eqs. (28), (48), (49), (34) and 
(43). The CDW C region is so narrow that it can hardly be 
distinguished on the total phase diagram. The dot line is a 
smooth link between the CDW0-CDW2C boundary near T c 
and its value H c = 2A /tt w 1.12T c0 at T = 0. 



CDW -CDW 2c boundary may be obtained qualitatively 
at arbitrary temperature as a smooth link between Fig. 
3 and the point H c = 2A /7r w 1.1T c0 (sec Fig. 5). 

In the case of non-uniform LOFF state the susceptibil- 
ity calculation in metal phase only gives the length of the 
optimal wave vector <7 op t| of the order-parameter mod- 
ulation. To determine the most energetically favorable 
combination of q : \q\ = \q op t\ one has to consider the 
third-order terms in |A| in the consistency equation. To 
determine which phase, CDW C or CDW2 C , wins we per- 
formed the similar procedure. Our procedure takes into 
account some peculiarities of the CDW state such as the 
influence of the " antinesting" harmonic in the electron 
dispersion. 

We have already mentioned the resemblance of the 
phase diagram obtained for CDW (Fig. 3) with that 
of layered superconductors with Rashba term. 41 The 
CDW0-CDW2C transition at low temperature also has 
many common features with the transition from uni- 
form superconducting state to LOFF state in layered 
superconductors. 42 In that case at the first critical field 
H c i the formation of a single soliton kink also becomes 
energetically favorable, and the functional shape of this 
soliton is well approximated by A(x) = Ao tanh(a;/a;o) 
both in CDW2C and LOFF states at low temperature. 
However, there are some important differences between 
these two systems. In the case of superconductivity 
the phase transition to the LOFF state was predicted 
to be of the second kind 42 , and the concentration of 
soliton walls increases gradually from zero with the in- 
crease of the magnetic field. In Ref. [17] the first- 
order phase transition was predicted between the uniform 
CDW and soliton phases. Our calculation also suggests 



the first-order transition between CDWo and CDW2 C 
phases. The first-order phase transition from metal to 
the LOFF phase was predicted by both the mean-field 42 
and renormalization-group studies. 7 The renormaliza- 
tion group analysis suggests 7 that the enhanced fluc- 
tuations of the LOFF state, associated with additional 
broken symmetries, are important to make this transi- 
tion of the first order. We do not consider the effect of 
fluctuations. According to the criterion in Refs. [5,6] the 
effects of fluctuations must not destroy the mean-field 
solution in organic metals where the interlayer transfer 
integral tb is, usually, much larger than T c q. The transi- 
tion from the normal state to the CDWj; phase observed 
on experiments 26-36 seems also to be of the second order. 
Further analysis of the CDWo-CDW x phase diagram in 
the case of imperfect nesting at other dispersion relation 
can give more accurate description and even new quali- 
tative features. 

Our results arc aimed to help to analyze the rather 
complicated phase diagram of CDW state in high mag- 
netic field observed in a number of organic metals, such 
as a-(BEDT-TTF) 2 KHg(SCN) 4 26 - 34 (see also a discus- 
sion in [38]), (Per) 2 M(mnt) 2 (M being Au,Pt,Cu) 35 ~ 37 
etc, where the transition from CDWo to high-field CDW^ 
phase is within the experimentally achievable magnetic 
fields. The above investigation may also be applied to an- 
alyze the phase diagram in nonorganic compounds with 
CDW ground state 2 (the attainable magnetic field in 
pulsed magnets is already comparable to the transition 
temperatures in some of these compounds) . An accurate 
quantitative description and comparison with the experi- 
mental observations in Refs. [26-36] may require the sub- 
stitution of a more realistic electron dispersion relation 
for each compound in the formulas derived above. With 
the two-harmonic tight-binding dispersion our study al- 
ready explains several qualitative features of the CDW 
state in high magnetic field. First, we show that the tran- 
sition from CDWo to CDWj state is of the first kind and 
is accompanied by the substantial jumps of the energy 
gap in electron spectrum and of the nesting vector that 
results in strong hysteresis. This fact being observed in 
many experiments on a-(BEDT-TTF) 2 KHg(SCN) 4 has 
not received a theoretical substantiation before. Second, 
we propose and prove the appearance of double-cosine 
CDW above the kink field which is consistent with the 
low temperature soliton solution. Third, we describe 
some properties of the cosine and double-cosine phases. 

To summarize, we have developed a mean-field theory 
of the CDW state in magnetic field which takes into ac- 
count imperfect nesting and the shift of the CDW wave 
vector. This allows a detailed study of the CDW prop- 
erties and the phase diagram in high magnetic field be- 
low the transition temperature. Our analysis gives the 
link between the previous theoretical results based on the 
susceptibility calculation 24 and the soliton solutions at 
zero temperature. 16 ' 17 Although our study of the CDW 2c 
state is only applicable close to the transition line T C (H) 
where the ration A a /T c can be considered as a small 
parameter, even an investigation in this region allows 
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to make some important conclusions about the struc- 
ture of CDW at high magnetic field. We have shown 
that the CDW X state at high magnetic field has predom- 
inantly a double-cosine modulation with wave vectors 
Q± = Qo ± q x . At perfect nesting the one-cosine CDW 
with shifted wave vector exists (according to the mean- 
field calculations) only in a very narrow region near the 
tricritical point (see Fig. 3 and Fig. 5). The finite sec- 
ond (" antinesting" ) harmonic in the electron dispersion 
(2) does not change this picture considerably. However, 



more specific dispersions or other perturbations may sub- 
stantially change the phase diagram. We also suggest 
an interesting peculiarity of the cosine phase — its spin 
asymmetry that allows to make a controllable spin cur- 
rent. 
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